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ABSTRACT 

The number density of rich galaxy clusters still provides the most robust way of normal- 
izing the power spectrum of dark matter perturbations on scales relevant to large-scale 
structure. We revisit this constraint in light of several recent developments: (1) the 
availability of well-defined samples of local clusters with relatively accurate X-ray tem- 
peratures; (2) new theoretical mass functions for dark matter haloes which provide a 
good fit to large numerical simulations; (3) more accurate mass-temperature relations 
from larger catalogues of hydrodynamical simulations; (4) the requirement to consider 
closed as well as open and flat cosmologies to obtain full multi-parameter likelihood 
constraints for CMB and SNe studies. We present a new sample of clusters drawn from 
the literature and use this sample to obtain improved results on erg, the normalization 
of the matter power spectrum on scales of 8 hr 1 Mpc, as a function of the matter den- 
sity and cosmological constant in a Universe with general curvature. We discuss our 
differences with previous work, and the remaining major sources of uncertainty. Final 
results on the normalization, approximately independent of power spectrum shape, can 
be expressed as constraints on <r at an appropriate cluster normalization scale i?ci- We 
provide fitting formulas for Rci and a(Rci) for general cosmologies, as well as for <Tg as 
a function of cosmology and shape parameter T. For flat models we find approximately 
08 — (0.495^Q'g37) f° r r = 0.23, where the error bar is dominated by uncertainty 

in the mass-temperature relation. 
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1 INTRODUCTION 



In theories of hierarchical structure formation the class of 
objects most recently formed holds a special significance. 
Observationally this class of objects is clusters of galaxies - 
the largest virialized structures in the present day universe. 
The local abundance of rich clusters of galaxies provides a 
strong constraint on the fluctuations in the matter density 
on scales of order 10 Mpc (Evrard 1989; Frenk et al. 1990 
Bond & Myers 1991; Henry & Arnaud 1991; Kaiser 1991 
Lilje 1992; Oukbir & Blanchard 1992; Bahcall & Cen 1993 
Hanami 1993; White, Efstathiou & Frenk 1993). Consistency 
with this constraint is one of the most important tests a 
model can pass, since the constraint is directly on the linear 
theory power spectrum, at a scale where there is an abun- 
dance of data. By fixing the normalization at wavelengths 
much smaller than those probed by COBE one obtains an 
accurate local normalization on scales relevant to much of 
structure formation, a long lever arm for constraining the 
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shape of the power spectrum, and a normalization to mat- 
ter fluctuations which is independent of galaxy bias. 

There have been several recent and detailed studies of 
the cluster abundance, including Bond & Myers (1996), Eke, 
Cole & Frenk (1996), Viana & Liddle (1996; 1999), Co- 
lafrancesco, Mazzotta & Vittorio (1997), Kitayama & Suto 

(1997) , Eke et al. (1998), Pen (1998), Wang & Steinhardt 

(1998) , Donahue & Voit (1999) and Henry (2000). However, 
even more recently, there have been several developments 
in terms of both theory and observation which suggest it 
would be useful to revisit this constraint. Firstly the addi- 
tion of ASCA temperatures (Tanaka, Inoue & Holt 1994) 
means that there is now a well defined local temperature 
function for clusters, with relatively small errors in temper- 
ature. Secondly a number of large N-body simulations have 
accurately determined the mass function of virialized haloes 
(e.g. Governato et al. 1999), finding non- negligible devia- 
tions from the old Press-Schechter (1974; hereafter PS) the- 
ory. For example the extremely large N-body simulations 
of the Virgo consortium have highlighted systematic depar- 



2 E. Pierpaoli, D. Scott & M. White 



tures from the PS predicted mass functions (Jenkins et al. 
2001), which alter the constraints on the power spectrum 
normalization coming from the cluster abundance. Thirdly, 
more ambitious hydrodynamical simulations of cluster for- 
mation (e.g. Frenk et al. 1999, and references therein) have 
resulted in improvements in the relationship between mass 
and temperature and a better estimate of its scatter. Finally, 
the increased sophistication of multi-parameter cosmologi- 
cal studies (in particular driven by recent CMB anisotropy 
measurements) requires that cosmological models with gen- 
eral curvature be considered (see e.g. White & Scott 1996). 

With these refinements our results are an improvement 
over other studies of the past few years. We point out ex- 
plicitly where we differ from other work, and also where we 
think things could be further improved in the future. The 
outline of the paper is as follows: we review some of the 
appropriate theory in §2; our local sample of clusters is pre- 
sented in §3; we describe our statistical method in §4; and 
we present our results and conclusions in §5 and §6. 



2 THEORY 

The abundance of rich clusters is tied to the normalization of 
the power spectrum (extrapolated to the present using lin- 
ear theory) through the Press-Schechter (1974) theory and 
its extensions (see Sheth, Mo & Tormen 2000 and references 
therein). The theory has always had a somewhat weak ana- 
lytic justification, its widespread adoption arising from the 
dual facts that it is easy to use and provides a remarkably 
good fit to more computationally expensive simulations. 

Although galaxy velocity dispersion and gravitational 
lensing mass estimates exist for many clusters (Carlberg et 
al. 1996; Smail et al. 1997; Girardi et al. 1998a; Allen 1998; 
Wu et al. 1998; Girardi et al. 2000), the direct estimation 
of mass through determination of the X-ray temperature is 
far more reliable - although the situation is certainly im- 
proving, allowing for estimates of the mass function (Gi- 
rardi et al. 1998b; Borgani et al. 1999). In the X-ray band 
there is considerably more high quality data on luminos- 
ity than on temperature, but there is enormous uncertainty 
in deriving mass M from luminosity Lx- Hence our point 
of comparison between theory and observation will be the 
temperature function of rich clusters, i.e. the number den- 
sity of clusters within a temperature range dT about T. This 
needs to be observationally determined over some range of T 
sufficiently high that gravitational physics dominates. Ob- 
servational and theoretical considerations place this limit 
at Tx ^ 3.5 keV (e.g. Finoguenov et al. 2000; Nevalainen, 
Markevitch & Forman 2000). Ideally we could compare the 
data with a temperature function estimated directly from a 
series of large cosmological hydrodynamic simulations (see 
e.g. Pen 1998), which included all of the physics relevant 
to determining the emission weighted inter-galactic medium 
(IGM) temperature of clusters. This however is currently 
computationally infeasible. Instead we shall make use of the 
fact that rich clusters are large virialized structures domi- 
nated by gravity, and thus we can factorize the problem and 
proceed in two steps. 

First we shall use the abundance of dark matter haloes 
of a given mass drawn from extremely large N-body simula- 
tions. Any dark matter halo large enough to host a cluster is 



unambiguously seen in such simulations. Next we shall use a 
mass-temperature relation, calibrated from hydrodynamical 
simulations, to convert from the (unobservable) virial mass 
to the IGM temperature, effectively using the larger volume 
of the N-body simulations to improve the statistics of the 
hydro simulations. Unfortunately the mass functions deter- 
mined from the N-body simulations are somewhat depen- 
dent on the method used to define haloes and their masses, 
and the precise definition of mass used there is not exactly 
what is used in the hydrodynamic simulations which cali- 
brate the M—T relation; the difference is expected to be 
small however (Jenkins et al. 2001). We discuss further de- 
tails of the M-T relation in §§2.4 and 2.5. 



2.1 The mass variance 

Our constraint will be on the variance of the density field, 
smoothed on some (comoving) scale R corresponding to a 
mass M = (47r/3)pi? 3 where p is the background density. In 
terms of the power spectrum 



a 2 (R,z) 



dk 

^ A 2 (k,z)W 2 (kR), 



(i) 



where A 2 = k 3 P(k, z)/(2tv 2 ), P(k) = \S k \ 2 is the matter 
power spectrum and W(kR) is the window function corre- 
sponding to the smoothing of the density field (see e.g. Pee- 
bles 1993). Our mass functions are fitted assuming a spher- 
ical top-hat smoothing, so 



W(kR) = 



Sji(kR) 
kR 



(2) 



where j\(x) is the spherical Bessel function of order 1. We 
are interested in both the normalization of the power spec- 
trum, for which we shall use as = <r(8 h~ x Mpc), and its 
shape. We use the Cold Dark Matter family of power spec- 
tra and parameterize the shape by F in the fitting formula 
of Bardeen et al. (1986). While the form of Eisenstein & 
Hu (1999) provides a slightly better fit to the shape, we will 
mostly quote the results in a F independent way, rendering 
this distinction unimportant. 

We write a(R,z) — a(R, 0) g(z)/g(0) where the growth 
factor g(z) can be computed numerically (Heath 1977; Car- 
roll, Press & Turner 1992; Cohn 1999; Hamilton 2001): 



, . 5 Qm da f a , , f da' 
9{z) =2 — d-r! da [dF 



(3) 



where the scale factor a = (1 + z) 1 , and the dimensionless 
time t = Hot. The Friedmann equation gives 



(4) 



where Ho = (a/a)t is the Hubble constant. The usual sym- 
bols Qm, = A/3-ffo and Ok = 1— Om — ^a are the density 
parameters (p/p C rit, with p cr it = 3Hq /8tvG) in matter, cos- 
mological constant and curvature, respectively. Photons and 
other relativistic species can be safely ignored. 
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Figure 1. A comparison of the Press-Schcchtcr form for the mul- 
tiplicity function, used in all previous work, with the Sheth & Tor- 
men (1999) model, which better fits the Virgo simulation results 
(see text). Note that Press-Schcchtcr systematically underesti- 
mates the abundance of the most massive cluster haloes, v ~ 10, 
and can either overestimate or underestimate the abundance of 
more typical clusters with v ~ few. 



2.2 Press- Schechter and modifications 

Within the PS theory and its extensions the (comoving) 
number density of objects of (virial) mass M is a function 
only of a{M). Using the scaled variable 



a{M) 



(5) 



with 8 C = 1.686, we write the mass function in terms of the 
'multiplicity function' f(v) as 

M 2 dn dlnM 



(6) 



p dM dlnv 

In the spirit of using PS as a fitting function to N-body 
simulations we do not imbue S c with a cosmology depen- 
dence (see next section) but rather keep it constant. The 
Press-Schechter model for f(v) is 



-i//2 



(7) 



and this formula has been extensively used to make pre- 
dictions of cluster abundances in different cosmologies. Re- 
cently Sheth & Tormen (1999) have proposed a correction to 
the Press-Schechter formula, motivated by a model of non- 
spherical collapse, that better fits large N-body simulations: 



vf{v) (x {l + v- p )v 1/z e 



1/2-0/2 



(8) 



where p = 0.3 and v — 0.707m. The normalization is fixed 
by the requirement that all of the mass be in haloes, i.e. 



/ 



f{y)du = 1 ; 



(9) 



for the above choice of parameters the normalization factor 
is 0.2162. Jenkins et al. (2001) have shown that Eq. (8) is in 
very good agreement with the simulations of the Virgo Con- 
sortium, except for very rare objects which shall not be of 



interest in this work. We compare the multiplicity functions 
in Fig. 1. Note that Press-Schechter systematically overes- 
timates the abundance of objects of cluster mass having 
I'&v'&Q (ff(.R) % 0.7) and underestimates the abundance 
of objects corresponding to a higher a(R). 



2.3 Spherical top-hats 

The spherical top-hat ansatz (Peebles 1993; Liddle & Lyth 
2000; Peacock 1999) models the formation of an object by 
the evolution of a spherical overdense region embedded in a 
homogeneous 'background' of mean density p. This region 
begins by expanding at the same rate as the background, but 
since it is positively curved the expansion slows, comes to a 
halt and the region collapses. Mathematically the evolution 
proceeds to a point of zero radius, however physically we 
assume that virialization occurs at twice* the turn-around 
time, resulting in a sphere of half the turn-around radius. 
The overdensity (relative to the background) at turn-around 
is 97r 2 /16 for an Einstein-de Sitter model. At virialization 
the background has become less dense and the sphere's den- 
sity grown by a further factor of 8 - we then denote the 
overdensity relative to the critical density by A c . This pa- 
rameter has the value 187T 2 in an Einstein-de Sitter model, 
and will in general be a function of Qm and f^A. The extrap- 
olation from linear theory of this overdensity is normally 
denoted S c . It is (3/20)(12tt) 2/3 ~ 1.686 for Einstein-de Sit- 
ter and varies by a few percent in other cosmologies. This 
density is used as a threshold in PS theory and its exten- 
sions. We shall neglect the small cosmology dependence and 
simply take 8 C fixed throughout. 

The value of the density contrast at collapse, A c , on 
the other hand, should be calculated for each model, since 
it will be important for the M—T relation. We computed A c 
by numerically integrating the equations of motion for the 
spherical top-hat collapse, including the correction to the 
virial theorem from the Ar potential (Lahav et al. 1991; 
§4.2). Our results can be fit to 2 per cent over the range 
0.2 < Q.m < 1.1 and < ft A < 1 by 



fin 



^2 C H 

i,j=0 



(10) 



where x = Qm — 0.2, y = f2 A and the coefficients dj are re- 
ported in Table 1. As an example we plot A c vs Qm for flat 
models, where the dashed line shows our fit and the solid 
line is the exact relation. We also checked this fit with the 
one provided by Eke, Navarro & Frenk. (1998b) and found 
an agreement at the 1 per cent level. Note that some authors 
use a different convention in which A c is specified relative to 
the background matter density - our A c is Qm times theirs. 



2.4 Halo mass definitions 

Before we turn to the mass-temperature relation we note a 
few important details about determining the mass of a dark- 
matter halo. Unfortunately there is no unique algorithmic 
definition of a dark matter halo, even within a 3D simulation 



* In the presence of a cosmological constant there is a small mod- 
ification to this. 
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Table 1. Coefficients Cij of the fitting formula, Eq. (10), for the 
collapse overdensity A c . 
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Figure 2. Showing our fit to the overdensity at virialization, 
A c vs !!m for flat models, where the dashed curve is our fitting 
function and the solid curve is the exact result. 

itself. Some particular group finders are in common use, but 
no single group finder is always used. For large objects such 
as clusters all group finders should be able to find all of the 
clusters, so this is not of immediate concern, although the 
degree of substructure will be highly variable between group 
finders. 

More disconcertingly there are a wide number of def- 
initions of halo mass in the literature, and they can differ 
by a large amount. Even for the mass-temperature relations 
calibrated by hydro simulations different authors use differ- 
ent definitions of 'mass', with the differences dependent on 
the cosmological parameters. In this sub-section we briefly 
review the relevant mass definitions and discuss which is the 
most appropriate for the mass function of Eq. (8). 

Although other halo finders are in common use, we 
shall deal exclusively with haloes found using the Friends-of- 
Friends (Davis et al. 1985) algorithm, hereafter called FOF. 
The FOF algorithm has one free parameter, b, the linking 
length in units of the mean inter-particle spacing. Commonly 
used values of 6 are 0.1, 0.15 and 0.2. The mass of the halo 
is simply the sum of the masses of the particles identified 
as part of the halo. An alternative (and more easily inter- 
preted) procedure is to use FOF to find candidate haloes, 
identify a halo centre (e.g. the centre of mass of the halo 



or, more robustly, the position of the most bound particle) 
and then to calculate the mass from the spherically averaged 
density profile about that centre. This is the technique typ- 
ically used to define the mass in hydrodynamic simulations 
which calibrate the M — T relation. 

In this spirit we define Ma as the mass contained within 
a radius rA, inside of which the mean interior density is A 
times the critical density: 

J r 2 dr p(r) = -p crit ri . (11) 

The 'virial mass' from the spherical top-hat collapse model 
would then be simply Ma c . Other masses in common use 
are M500 and M200 where the latter is approximately the 
virial mass if Qm = 1. 

For large mass haloes many of these definitions are re- 
lated on average simply by a factor. We can estimate this 
factor by assuming that haloes have a universal profile, for 
example the NFW form (Navarro, Frenk & White 1996): 

p(r) ocaT^l + x) -2 , (12) 

where x = r/r s and r s is a scale radius usually specified in 
terms of the concentration parameter c = r2oo/r s - Navarro 
et al. (1996) refer to r2oo and M200 throughout as the 'virial 
radius' and 'virial mass' respectively. Again, N-body simula- 
tions have shown that the concentration parameter is a weak 
function of virial mass, having the value c ~ 5 for masses 
characteristic of clusters. We can use this profile to relate 
the various mass definitions, as shown in Fig. 3. 

Unfortunately, while this model works well for convert- 
ing between spherically averaged mass definitions based on 
Ma , there is a large scatter for masses based on group mem- 
bership, such as Mfofo.2- This is because with b — 0.2, FOF 
can link together neighbouring haloes in supercluster-like 
structures, increasing the mass assigned to the structure 
compared to the Ma estimators. Thus the 'FOF' lines in 
Fig. 3 must be taken as highly uncertain. Comparing the 
mass function from a high-resolution N-body simulation (of 
the Ostriker & Steinhardt 1995 concordance model) with 
the universal form of Jenkins et al. (2001) we find the best 
match is obtained if we interpret their mass as Ma c . We 
shall assume below that this remains true independent of 
cosmology. But we note that this ambiguity in the defini- 
tion of mass remains a significant source of uncertainty. 

2.5 The M T relation 

The mass function of rich clusters is itself not observable. 
However the local temperature function is reasonably well 
known (see §3). To predict the latter from the former we 
need a relation between emission weighted IGM temperature 
and cluster virial mass. Our results are quite sensitive to 
the choice of M—T relation, and currently the uncertainty 
in this relation is the largest theoretical source of error in 
determining as (see also Voit 2000). 

Recent observational determinations of the mass- 
temperature relation (Horner, Mushotsky & Scharf 1999; 
Nevalainen, Markevitch & Forman 2000) disagree at the sev- 
eral tens of percent level (in mass at fixed temperature) 
when using different estimators of the cluster mass. The 
virial mass of a cluster is a notoriously difficult quantity to 
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Figure 3. Relations between various definitions of the mass of a 
halo as a function of fin assuming the halo density profile follows 
the NFW form with concentration parameter c = 5. The FOF 
based masses are only crude approximations in this model (see 
text). 



Table 2. The mass-temperature relation determined from hydro- 
dynamical simulations. The quoted value of f3 is that relevant for 
use in Eq. (13). From top to bottom the references are: EMN 
(Evrard ct al. 1996); ENF (Eke et al. 1998); BN (Bryan & Nor- 
man 1998); YJS (Yoshikawa, Jing & Suto 2000); TC (Tittlcy & 
Couchman 2000) and Tctal (Thomas ct al. 2000). We have quoted 
emission weighted temperatures where available, while the last 
two values (marked with an asterisk) are core temperatures. 



Name 








h 


13 


EMN 


1.0 


0.0 


0.10 


0.50 


1.21 


EMN 


0.2 




0.10 


0.50 


1.42 


ENF 


0.3 


0.7 


0.04 


0.70 


1.33 


BN 


1.0 


0.0 


0.06 


0.50 


1.10 


BN 


1.0 


0.0 


0.10 


0.65 


1.04 


BN 


1.0 


0.0 


0.08 


0.50 


1.04 


BN 


0.4 


0.0 


0.06 


0.65 


1.08 


YJS 


0.3 


0.7 


0.03 


0.70 


1.48 


TC* 


1.0 


0.0 


0.10 


0.65 


1.61 


Tetal* 


1.0 


0.0 


0.06 




1.23 



obtain observationally with high accuracy - while different 
estimators clearly correlate well, they disagree at the level of 
accuracy required here. Such observations do however pro- 
vide general support for the functional form and scalings 
predicted by the spherical collapse model (see below) for 
clusters of sufficiently large mass. Similar scalings are seen 
in hydrodynamic simulations of galaxy clusters in a cosmo- 
logical context. In the simulations the total mass of a cluster 
is easy to obtain (though convention dependent, §2.4) and in 
what follows we shall use a M—T relation derived from sim- 
ulations. While these hydrodynamic simulations show good 
agreement for the total mass and X-ray temperature prop- 
erties of clusters (Frenk et al. 1999) there are several un- 
certainties which enter when comparing the simulations to 
observations and which are important to note. 

As with all simulation derived results there are issues 
related to numerical convergence. With the latest round of 
high resolution simulations the situation in this respect has 
improved dramatically. However, the spectrally measured 
cluster temperatures may not coincide with the mass or 
emission weighted temperature estimated from the simu- 
lations, due to the influence of soft line emission (Math- 
iesen & Evrard 2001). Secondly, most simulations are done 
with purely adiabatic hydrodynamics, which ceases to be 
a good model for the lower mass/temperature clusters. 
It is also possible that energy injection (possibly from 
SNe, AGN or galaxies) has altered the M—T relation, 
again an effect thought to operate preferentially on lower 
mass/temperature clusters. Thirdly the simulations usually 
predict the emission weighted temperature, which only con- 
verges if data are used out to a radius ~r 5 oo- Observers 
often probe different ranges of radius in determining the 
cluster temperature and perform more sophisticated mod- 
elling, which can lead to discrepancies between the measured 
and simulated temperatures. Finally, there still remain sig- 
nificant (for our purposes) calibration uncertainties for the 
detectors. 



With these caveats in mind, the results from the simu- 
lations can be quoted in terms of corrections to the spherical 
collapse model which relates the mass to the (virial) temper- 
ature of the hot IGM. For an object virialized at a redshift 
z we have^ 



M(T,z) 
lQiB/i-iMo 



1-2 



«a(z) 



-3/2 



(13) 



where T is in keV, A c is the mean overdensity inside the 
virial radius in units of the critical density and, from Eq. (4), 
E 2 = Q M (1 + z) 3 + Oa + Ok(l + zf. Note that A c is a red- 
shift dependent variable, and should be evaluated using the 
appropriate CIa(z) and fl M {z). The term in square brackets 
is a correction to the virial relation arising from the addi- 
tional r 2 potential in the presence of A (Lahav et al. 1991; 
Viana & Liddle 1996; Wang & Steinhardt 1998). It provides 
only a small correction and, though we include it, it can be 
neglected at the present level of accuracy. 

Ideally the normalization and scatter of the M—T rela- 
tion would be determined by simulations, so that theoretical 
models can be compared with the data rather more directly 
(e.g. Pen 1998). However Eq. (13) is a remarkably good fit 
to the simulations, which are sufficiently computationally 
demanding that they cannot explore parameter space effi- 
ciently. Thus we rely on a hybrid approach where the coef- 
ficients are determined from simulations, while the scalings 
are taken from simple theoretical models (Mathiesen 2000). 
Specifically we use the hydrodynamic simulations to deter- 
mine [3, together with Eq. (13) for the mass and redshift 
dependence. In practice we used the M—T relation at an 



t Our definition of /3 differs from that of Henry (2000): /^Hcnry = 
1.42//?, and our f3 should not be confused with the slope of the 
emissivity profile of clusters. 
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observed z = 0.053, corresponding to the median redshift 
of our cluster sample (see §3). We then shift the resulting 
a(z = 0.053) value to z = using the growth rate, Eq. (3), 
which is a significant 5 per cent correction. 

The simulations normalize j3 assuming that the virial- 
ization redshift is the redshift of observation. In principle one 
could attempt to correct for the virialization redshift depen- 
dence, however we have chosen not to do this, and here we 
differ from Viana & Liddle (1999) and Wang & Steinhardt 
(1998), for example. The simulations, which clearly include 
the full effects of variations in the virialization redshift, give 
an M—T relation well fit by Eq. (13) if z is interpreted as 
the redshift of observation. We believe that the effect of dif- 
fering virialization redshifts is included in the simulations 
as part of the scatter about the mean relation (see below) , 
and so to add an additional effect by hand would be incor- 
rect. Comparison of the scatter in the M—T relation found 
by Bryan & Norman (1998) in a full cosmological simula- 
tion with that of Evrard, Metzler & Navarro (1996), who 
use constrained realizations (and thus constrained forma- 
tion times), suggests that in fact the effect of scatter in the 
virialization redshift is a very small source of the total scat- 
ter in the relation. Most of the scatter is due to the different 
merging histories (see also Cavaliere, Menci & Tozzi 1999). 
While further simulations will be needed to address this is- 
sue properly, recent work (Mathiesen 2000) suggests that 
minor merger events have a more important influence on 
the evolution of the temperatures than major mergers (and 
hence formation time). 

A summary of recent numerical experiments which con- 
strain (3 is given in Henry (2000). We have taken Henry's list, 
added some recent work and corrected one of the relations 
to use M v i r . Our results are shown in Table 2. The M—T 
relation of Evrard et al. (1996) defines mass as M200, so that 
in the context of the spherical top-hat model there is pre- 
dicted to be an £2m dependence to the prefactor of the scaling 
relation. Correcting for this scaling brings the /3s obtained 
in these simulations into better agreement, but they still 
disagree slightly, suggesting that the scaling is only approx- 
imately observed. We postulate that this is because Qb/£Im 
changes drastically between the simulations. The Ms of Eke 
et al. (1998), Bryan & Norman (1998), and Yoshikawa et 
al. (2000), on the other hand, are already the 'virial' mass 
in the sense of the spherical top-hat model. Also quoted in 
Table 2 are the /3 values from Tittley & Couchman (2000) 
and Thomas et al. (2000). While all of the other tempera- 
tures in the Table are emission weighted, these authors use 
the average temperature within a core region, which could 
be slightly different. 

As shown in Table 2 the values of (3 spread from near 1 
up to 1.6 with no obvious peak of preferred values. We adopt 
the mean opinion on this issue, allowing for (3 = 1.3, with a 
10 per cent systematic error (i.e. variation between the simu- 
lations) in the mass, plus a 10 per cent statistical error. The 
second error arises from the intrinsic scatter in individually 
determined M—T relations, and is mostly due to the merg- 
ing history of the clusters. Most of the simulations agree on 
the scatter about the mean relation quite well, though Eke 
et al. (1998) find a slightly enhanced scatter compared to 
the other authors. This dispersion is accurately modelled as 
a Gaussian, and we treat the systematic uncertainty as a 
Gaussian also. The /3 value advocated by Henry (2000) is 



slightly lower, (3 = 1.17, with a suggested systematic error 
of 4.1 per cent, while the value we would obtain using the 
observational determinations is close to unity. 

2.6 Summary of modelling 

In summary, we use the mass function of Jenkins et 
al. (2001), interpretting the mass as the top-hat virial mass. 
For an object at fixed mass we randomly assign a temper- 
ature using Eq. (13) where A c is given by Eq. (10) and 
the M—T relation is evaluated at a median redshift inter- 
preted as the redshift of observation. From Table 2 we choose 
j3 — 1.3±0.13±0.13 with Gaussian errors. Here the first error 
represents the scatter about the mean relation and the sec- 
ond is the systematic uncertainty in the prefactor from the 
different calculations. We carry out the comparison between 
theory and data at the median redshift of the observations, 
and then correct the final normalization to the appropriate 
z — value. 



3 DATA 

3.1 Definition of cluster sample 

There has been progress recently in observationally deter- 
mining the temperature function of nearby clusters, and we 
have compiled from the literature a new sample with which 
to constrain the normalization of the power spectrum at 
z~0. 

The sample is adapted from the 30 clusters compiled 
by Markevitch (1998). His sample was selected from bright 
clusters having flux above 2 x 10~ 14 Wm~ 2 in the ROSAT 
0.2-2 keV band, and over the redshift range z — 0.04-0.09. 
The flux limit is a factor of 4 above the nominal flux limit of 
the ROSAT Brightest Cluster Sample (Ebeling et al. 1998), 
and hence the sample is expected to be close to complete 
for these fluxes. The upper redshift limit is imposed by the 
lack of clusters bright enough to be detected, and the lower 
limit by ROSAT selection effects. Clusters are also excluded 
with galactic latitude |6| < 20°, where observations become 
affected by the Galaxy. This sample is close to volume lim- 
ited at the high-temperature end, and the numbers can be 
corrected for incompleteness at the low-temperature end by 
using the effective volume (see §4.1). Markevitch (1998) de- 
termined temperatures by excising the central regions from 
clusters to approximately correct for the effects of cooling 
flows, and used a hybrid approach of combining ROSAT 
emissivity profiles with ASCA data to fit Tx. As we discuss 
below, we believe that somewhat better temperatures are 
now available based on careful fitting of ASCA data alone. 
Using mainly these temperatures, together with a somewhat 
different selection, we have attempted to define an effectively 
temperature selected sample. 

White & Buote (2000) use a maximum likelihood 
method to determine radial temperature profiles for clus- 
ters using ASCA data alone. This includes a Monte Carlo 
method for redistribution of X-ray photons by ASCA's com- 
plex optics. They account for cooling flows using a single 
extra parameter in their fits, and find that the cooling flow 
corrected temperatures are generally consistent with what 
would be fitted to the outer regions of the clusters, although 
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with less uncertainty. They also do not find the temper- 
ature fall-off at large radii found by Markevitch (1998), 
about which there has been much discussion in the literature 
(e.g. Irwin, Bregman & Evrard 1999). We regard the temper- 
atures presented in White (2000) as representing the most 
careful analysis of X-ray temperatures from available ASCA 
data. Independent BeppoSAX data, available for about a 
quarter of these clusters give temperatures in good agree- 
ment (Irwin & Bregman 2000). These values are unlikely to 
improve significantly until Chandra and XMM-Newton data 
become widely available. 

We have constructed our cluster sample in a similar way 
to that given by Markevitch (1998), although we use tem- 
peratures from White (2000) when available, supplemented 
by other temperatures from the literature. The sample of 
course has a large overlap with other low-redshift cluster 
samples, such as those of Edge et al. (1990), Henry & Ar- 
naud (1991), Henry (2000) and Blanchard et al. (2000) - 
although we neglect the lowest redshift clusters for reasons 
of incompleteness, as well as concerns about biases intro- 
duced by sample variance. 

Since many of the White (2000) temperatures show 
significant changes compared with Markevitch (1998), we 
also need to reconsider the completeness of the sample. We 
use the estimated Lx — Tx relationship derived from Marke- 
vitch (1998) to decide whether a cluster would be above the 
ROSAT flux limit based on the improved value of Tx- Since 
we find that the White (2000) temperatures are a factor 
~ 1.14 higher than the Markevitch (1998) temperatures on 
average for the clusters in common, we correct the Lx—Tx 
relation by this factor. Explicitly we use 

L X = 1.07 X 10 37 (^) 2 V 2 W, (14) 

where as usual h = ffo/(100 kms _1 Mpc _1 ). Marke- 
vitch (1998) finds a relatively small scatter about this re- 
lationship once cooling-flow effects have been corrected for. 
(A steeper temperature dependence is found for bolometric 
X-ray luminosity or luminosity over higher energy ranges.) 
The use of this relationship allows us to avoid using any 
specific flux information for a particular cluster, which is 
advantageous since this information is somewhat uncertain, 
varying significantly between instruments and between anal- 
ysis methods. Provided Eq. (14) is approximately correct, 
the details will be a higher order correction to our er 8 con- 
straints. 

Our sample is thus effectively temperature-selected at 
each redshift, and hence we will only need to correct for 
the volume sampled at each temperature over the full red- 
shift interval. We use CMB-frame redshifts from the recent 
compilation of Struble & Rood (1999) for the Abell clus- 
ters, and the redshifts given in White (2000) or Ebeling et 
al. (1998) otherwise. At these low redshifts only a small er- 
ror is made by assuming non-expanding Euclidean space, 
since the correction to the luminosity distance is 0(z/4), 
which is ~ 2 per cent at worst. It is easy to compute the 
relationship exactly, although this would in principle have 
to be re-calculated for each model. We use redshift as an 
exact distance indicator, which is another reason to cut off 
the sample at low redshifts, where this will cease to be a 
good approximation. The effective volume, as a function of 
temperature, is shown in Fig. 4. Because of the large com- 




kT (keV) 

Figure 4. The effective volume of the sample as a function of 
temperature. Note that the sample is volume limited for clusters 
above about 7.5 keV, and that the correction factor is not too 
drastic for kT Jt 3.5 keV. 



pleteness correction required at low temperature, we make 
a further cut of all clusters with Tx < 3.5 keV. This also 
corresponds approximately to a restriction to clusters with 
temperatures which are dominated by gravitational physics 
(Mohr & Evrard 1997; Balogh et al. 1999; Xue & Wu 2000). 
Our final sample of 38 clusters is presented in Table 3. 



3.2 Supplementary sample 

Since we will be considering the temperature errors in our 
Monte Carlo error analysis (§4.2), some clusters will scatter 
out of the selection cuts and hence it is important to include 
clusters which could scatter into the sample. We include in 
Table 4 additional clusters with Tx + 2a being above our 
3.5 keV cut-off, or above the temperature-derived flux limit 
at their redshift. This limit corresponds to 

TV / z \ °- 95 

keV> 2 - 65 (^) ' < 15 > 

and is indicated in Fig. 5 by the roughly diagonal line. To 
populate this supplementary sample (as well as to check 
completeness of the main sample) we extended the ROSAT 
selection down a further factor of 2 to 1.0 x 10~ 14 Wm~ 2 . 
In particular we scrutinized the BCS sample of Ebeling et 
al. (1998), together with the RASS1 Bright Cluster Survey 
of de Grandi et al. (1999) for the southern extension, and 
the overlapping XBACS sample (Ebeling et al. 1996) came 
from either White (2000) or the David et al. (1993) compi- 
lation when available. The REFLEX survey (see Bohringer 
et al. 2000) may ultimately be better for this purpose, but is 
not yet published. For some clusters we used temperatures 
from the White, Jones & Forman (1997) deprojection mod- 
elling of Einstein data. For a few cases we had to resort to 
estimates (Ebeling et al. 1998) based on the X-ray luminos- 
ity - for those cases we adopted a representative error of 
±1.00 keV, which flags them in Table 4. 
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3.3 Comparison with other samples 

Our final sample differs from that of Markevitch (1998) in 
several details. We chose to widen the sample down to z = 
0.03, since that increased the statistics, while we could find 
no evidence that there was any significant bias introduced. 
This can be seen in Fig. 5, where we have shown the clusters 
(solid points) together with our selection cuts in redshift and 
in temperature. Clusters which could scatter into our sample 
are indicated by open symbols. 

The White (2000) temperatures were generally a little 
higher than those of Markevitch (1998), resulting in several 
clusters entering our sample, explicitly Abell 193, Abell376, 
Abelll775, Abell 2255, Abell 3532 and IIZwl08. On the 
other hand Abell 780, Abell 1650 and MKW3s were too cool 
at their redshifts, while Abell 3112 and 2A0336 even failed 
the 2cr 3.5 keV cut. Abell 2199, Abell 2634 and Abell 4038, 
which are in some other samples, were lost here because 
their CMB-frame redshifts (Struble & Rood 1999) - as op- 
posed to the heliocentric redshifts more commonly given - 
are slightly less than 0.03, while Abell 3921 is higher than 
0.09. Furthermore we added Abell 496, Abell 576, Abell 2063, 
Abell 2107, Abell 2147, Abell 2151a, at 0.03 < z < 0.04. 

Our final sample of 38 clusters is presented in Ta- 
ble 3, where we list the common name, redshift (taken 
from Struble & Rood 1999 or White 2000), X-ray tempera- 
ture and la error bar. These temperatures are mostly from 
White (2000), but in a few cases from Markevitch (1998) 
errors scaled from 90 per cent confidence, or from the 80 
per cent confidence region estimates of White et al. (1997). 
Table 4 lists the additional clusters which could scatter into 
our sample within their temperature uncertainties. Through 
this selection process we believe we have a reasonably con- 
structed temperature-selected sample. 

Though we do not use it directly in the analysis (see 
§4.1) we have constructed the temperature function of our 
sample for comparison with earlier work. We show this in 
Fig. 6, where we have only used the clusters in Table 3. 
Some of the other estimates plotted were derived for specific 
cosmological models or at other redshifts, so they cannot be 
compared in great detail. However, it is clear that we are in 
general agreement with the temperature function estimated 
by Markevitch (1998), although a little higher. More striking 
is that both of these estimates are considerably higher than 
those of Henry (2000) and of Eke et al. (1996) and Viana 
& Liddle (1999) which are based on the earlier Henry & 
Arnaud (1991) sample. The main reason for this difference 
is the correction for cooling flows. We believe that the most 
physically realistic comparison of the z — clusters with the 
results of simulations is to correct for cooling flows, since the 
full cooling-flow physics is not contained in the simulations. t 
This gives rise to an added complication when carrying out 
evolutionary studies of the high-z vs low-z samples, since 
detailed information about cooling flows in high- z clusters 
tends not to exist. However, for our purposes it seems clear 
that we should use cluster temperatures which have been 



Table 3. Our nearby cluster sample. The data are adapted from 
the sample of Markevitch (1998) with White (2000) temperatures, 
and Struble & Rood (1999) redshifts, as discussed in the text. 
Reported errors are Icr. Note that the White (2000) method may 
overestimate the temperature for some clusters for which the ex- 
istence of a central cooling flow is debatable, e.g. A399, A401, 
A754, A1775, A3158 and A3562 (M. Markevitch, private com- 
munication). Removing these clusters has negligible effect on our 
results, and we prefer to use the White (2000) temperatures, when 
available, for the sake of consistency. 



Name z T x 



A85 


0.0543 


D -' 4 -0.50 


A119 


0.0430 


ft nc:+0.55 
D.UO_ 43 


A193 


0.0476 


4 20 +1 '°° 

^• zu -0.50 


A376 


0.0472 


5.70l?'?7 


A399 
A401 


0.0712 
0.0725 


O cc+l-92 
y - dd -0.96 

io.68±S;ii 


A478 


0.0869 


7 49+ ' 71 

' -^ z -0.54 


A496 


0.0317 


4 en +0.17 
4 - ol -0.15 


A576 


0.0377 


4 02+ ' 20 


A754 


0.0530 


12.85li"S 


A1644 


0.0461 


4 7n +0.49 


A1651 


0.0832 


7 n r+0.84 
' ■ lo -0.62 


A 1736 


0.0446 


4.02tg;£ 


A1775 


0.0705 


8.70t° 2 ;« 


A 1795 


0.0619 


7 26+ ' 51 

'• zo -0.40 


A2029 


0.0761 


o 9 n+0.58 
°- zz -0.20 


A2063 


0.0341 


o qn+ ' 51 


A2065 


0.0714 


6 1Q+ ' 70 
u - la -0.70 


A2107 


0.0399 


4 Q1+0.57 
^• di -0.35 


A 91 49 


D8Q7 


1 n qft+ 2 ' 56 


A2147 


0.0338 


. 4t;+o.5i 

°-^°-0.38 


A2151a 


0.0354 


3.80tg;™ 


A2255 


0.0794 


7.76li;° 01 


A2256 


0.0569 


8.69tl:g. 


A2589 


0.0402 


3.70li;° 04 


A2657 
A3158 


0.0390 
0.0585 


3.891°;- 

o qq+1.43 


A3266 


0.0577 


q 6 o+0 97 


A3376 


0.0444 


4.38l°;?« 


A3391 


0.0502 


6-901^ 


A3395 


0.0494 


4.801°;- 


A3532 


0.0542 


4 70+0- 39 


A3558 
A3562 
A3571 


0.0468 
0.0478 
0.0379 


6 60+ ' 50 
ft QR+ 1 ' 77 

u - yo -0.95 
o i 9+0.42 
°- 1 -0.39 


A3667 


0.0544 


8 11+0.82 
°- 11 -0.73 


A4059 


0.0463 


4 05+ ' 23 

^•LM_0.19 


IIZwl08 


0.0493 


4 40+ 100 
^•^ u -i.oo 



t The cooling flow correction applied by White (2000) may over- 
estimate the temperature for some clusters like A754. However, 
we decided to stick to the White (2000) results in all cases, for 
reasons of consistency. 
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Table 4. Additional clusters which could scatter above our se- 
lection criteria within their ±2<r temperature uncertainties. Ad- 
ditional temperatures come from David (1993) and White et 
al. (1997). Those with errors given as ±1.00 kcV are estimates 
from X-ray luminosities. 



Name 


z 




A133 


0.0554 


4 00 +14 ° 


A168 


0.0438 


2-60l i» 


A780 


0.0527 


4.49t°« 


A1650 


0.0833 


D.OO_ 43 


A1800 


0.0743 


O-J-U-l .00 


A1831 


0.0603 


4 on+ 100 


A2052 


0.0338 


3.30l°;i« 


A2061 


0.0772 


5-60±i;°g 


A2151 


0.0354 


3.04+g;« 


A2249 


0.0804 


5.60li;°° 


A2495 


0.0763 


5.00±i;gg 


A2572a 


0.0391 


2-80l°;« 


A2593 


0.0401 


q 1 n+i-50 

°- lu -0.90 


A2665 


0.0544 


4 onl" 1 ' 00 
^• zu -i.oo 


A2734 


0.0613 


A co+1-00 
*- ou -1.00 


A3528a 


0.0516 


4 in+ 1 -00 
^• lu -1.00 


A3560 


0.0477 


3 40 +100 


A3695 


0.0882 


°- lu -1.00 


A3716 


0.0450 


O 1Q + 1.00 

°' lu -1.00 


A3809 


0.0608 


4 •?n +100 


A3822 


0.0747 


5 50+ 1 ' 00 

J,JU -1 .00 


A3822 


0.0747 


cc rn+100 
O.OU_j oo 


A3880 


0.0572 


4 00+ 1 ' 00 


EXO0422 


0.0397 


2 Q0+ ' 50 


MKW3s 


0.0434 


o 71+O.I6 
°- ' -0.19 


RXJ1733 


0.0330 


2 60+ 1 00 

Z - DU _1.00 


S0405 


0.0613 


4 30+ 1 ' 00 
^• ou -i.oo 


S1101 


0.0580 


n m +1.20 


SC1327 


0.0476 


n q 9 +0.44 
°- az -0.32 


Zw235 


0.0830 


t- oo+l.OO 
°- zu _1.00 


Zw5029 


0.0750 


O.OU_j oo 


Zw8276 


0.0757 


r 7n + 1.00 


Zw8852 


0.0400 


q 07+O.IO 
- ' -0.10 



corrected as carefully as possible for the effects of cooling 
flows. 



4 STATISTICAL APPROACH 
4.1 The likelihood function 

Rather than fitting to a measure of n(> fcTx) at some fidu- 
cial Tx, we try to fit to all the data in a way which accounts 
for the Poisson distribution of clusters at each temperature, 
using a likelihood method which takes full account of the 
errors. Our approach is similar to that of some other studies 
(e.g. Eke et al. 1996; Markevitch 1998; Henry 2000; Blan- 
chard et al. 2000), but we describe it here in full, so that the 
differences can be appreciated. 

Given a set of data with which to compare (§3) we com- 
pute the likelihood function for any given theory as follows. 




Figure 5. An indication of our selection function (dashed lines). 
We chose clusters (solid symbols) with CMB-frame rcdshifts 
within 0.03 < z < 0.09 and a tempcraturc-dcpcndcnt function to 
account for a flux limit. We also cut out all clusters below 
fcTx = 3.5keV. However, we kept clusters (shown as open circles) 
which could scatter into our acceptance region within their ±2c 
temperature uncertainties (which arc plotted). Rcdshift errors are 
here assumed to be negligible. The lack of clusters at z ~ 0.065 
is a well-known effect caused by a genuine lack of supcrclustcrs 
at those distances. 




4 5 6 7 8 9 10 

kT x (keV) 



Figure 6. The temperature function constructed from the sam- 
ple in Table 3 (solid line), compared to some earlier estimates. 
The dashed line is from Markevitch (1998), the boxes from 
Henry (2000), the filled circles from Eke et al. (1998), and the 
open circle from Viana & Liddle (1999). The curved line is our 
best fit for a flat cosmology with f!|n = 0.3 
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We break the range of temperatures under consideration 
into a large number of bins, chosen to be narrow enough 
so that the probability of two clusters occupying the same 
bin is very small. Then for a given Monte-Carlo realization 
(§4.2) of the temperatures of the set of clusters we place the 
clusters into the appropriate temperature bins, giving an oc- 
cupation number r\i — or 1 for each bin. For a predicted 
number density Mi = n(T)dT, depending on our cosmolog- 
ical parameters, the mean occupation number of each bin 
is Hi = AfiVi <C 1. Here Vi is the volume of space to which 
clusters in bin i can be seen: 

Vi = ^(dl-d\), (16) 

if the survey solid angle is Q. s , which is 4tt(1 — cos 70°) ~ 
8.27 sr here. With a Euclidean assumption di — czi/Hq, 
while 

where /n m is the limiting flux of the survey (2 x 10 14 W m~ 2 ) 
and L(Tx) is obtained from the luminosity-temperature re- 
lation, Eq. (14), given the cluster's temperature. Then the 
likelihood of observing this combination of clusters is simply 

In £ = [(»?< - + r?i In (1 - cxp(-^))] , (18) 

i 

where by assumption only one of the two terms is non-zero 
for each bin i. This correctly accounts for the Poisson errors 
and uses the full temperature information from the sample. 

4.2 Monte Carlos 

Given that there are non-negligible uncertainties occurring 
in several places, that the calculation is quite non-linear, 
and that n(T) is a steeply falling function, it is important to 
treat errors carefully. The only reasonable way to do this is 
through a Monte Carlo approach (also emphasized by Viana 
& Liddle 1999 and Blanchard et al. 2000), which we now de- 
scribe. Firstly, we choose a temperature for each cluster (in 
Tables 3 and 4) by generating Gaussian random numbers, 
using the upper and lower error bars each 50 per cent of 
the time. Then we form a new sample by culling all clusters 
with temperatures which fail our selection cuts (as described 
in §3). Next we adopt a specific M—T relation drawn from 
the central value and systematic range discussed in §2.5, 
together with an additional scatter, different for each tem- 
perature considered, arising from the intrinsic scatter in the 
M—T relation. For each temperature, we use the M—T re- 
lation to calculate n(T) according to the mass function in 
Eq. (8). 

For each cosmological model considered, we then maxi- 
mize the likelihood to find a best-fitting normalization from 
which to determine as- This whole process is repeated 1000 
times to obtain a distribution of as values for each cosmol- 
ogy- 



5 RESULTS 

Our main result is the power spectrum normalization as a 
function of cosmological model. We quote the normalization 



Table 5. Coefficients of the fitting formula, Eq. (20), for the mean 
and errors of a(R C i) = p Q^+™ n ™+™ n ^ . 







parameter 






Po 


Pi 


P2 


P3 


mean 


0.505 


0.150 


-0.233 


0.048 


-la 


0.473 


0.124 


-0.234 


0.053 


+la 


0.532 


0.171 


-0.222 


0.044 



for the variance on a scale Rci that corresponds to a clus- 
ter of 6.5 keV forming now. This is advantageous because 
the value of a(Rc\) is almost independent of the power- 
spectrum shape F (see also Blanchard et al. 2000), while at 
fixed <r(i?ci) the value of as can vary by as much as 15 per 
cent as T spans the 68 per cent confidence range, 0.19-0.37, 
of Eisenstein & Zaldarriaga (2001). Fitting formulae accu- 
rate at the 1 per cent level for Rci and a(Rci) covering the 
ranges 0.2 < S7m < 0.8 and 0.3 < Qa < 1 are given by: 

Rci = Po « M (P1+P2 nM+P3 " a) /r'Mpc ; (19) 

a C i = a(R C i) = Po fi M (P1+P2 nM+P3 " A) • (20) 

Here p = (9.086, 0.354, 0.058, 0.049) gives the coefficients in 
Eq. (19), while those for Eq. (20) can be found in Table 5, 
along with fits to the ±1<t limits on a(Rci)- The validity of 
Eq. (20) has been tested also in the case of open (i.e. Qa = 0) 
and Einstein-de Sitter (i.e. Qm = 1) models, and for this 
wider range of parameters the maximum deviation from the 
fit is still only 2 per cent. For the concordance model (Om = 
0.3) we find Rci — 14.8 /i _1 Mpc. A slightly less precise, but 
simpler fit for flat models is given by 

Rci =± 9.0 nM°' 41 ft _1 Mpc. (21) 

This fits to better than 1 per cent over 0.2 < Qm < 1. 

For explicit constraints on as we verified that the fol- 
lowing equation is a good fit (with a maximum error of 2 
per cent), in the range 0.19 < T < 0.37 (Eisenstein & Zal- 
darriaga 2001), 0.2 < n M < 0.8 and 0.3 < fl A < 1: 

as = aci go ^a 1 Q M fe+ « nM) r (94+95nM) . (22) 

Here the vector of coefficients is given by q — 
(1.049,0.029,0.546,0.350,0.278,-0.312). For flat models 
(i.e. Qa = 1 — £^m) this corresponds approximately to 

as^ (0.495^7) n M °' 60 (23) 
for the value F = 0.23 explicitly, and for general F: 

a C i = 0.462_ 029 S2 M , (24) 

where the exponents in brackets refer to the best power- 
laws for the +la and —la limits. It is interesting to notice 
that, despite the many differences in our approach, this re- 
sult does not differ dramatically from other studies, such as 
those reported in Fig. 7. This is true even for those stud- 
ies with apparently quite discrepant temperature functions, 
such as Henry (2000) . As we discuss below this is partly due 
to the chance cancelation of several small changes affecting 
the derived as . 
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Figure 7. The normalization erg for flat cosmologies with shape 
parameter fixed at T = 0.23 as derived here (solid, red ) with 1— 
a limits (dashed red lines), compared with the fits from Viana 
& Liddle (1999) (dot-dashed dark blue line), Wang & Stein- 
hardt(1998) (blue 3 dots-dashed), Eke et al. (1996) (dotted 
green), and Borgani et al.(1999) (cyan). 



We performed several checks on our results. First we 
made sure that analysing the sample using a single (me- 
dian) redshift did not introduce a significant bias. We split 
the sample into low-z and high-z halves and analysed them 
independently. The variation in aci between the low-z, high- 
z and combined samples was only 2 per cent for a flat model 
with Qm = 0.3. We also gradually eliminated clusters with 
redshift below a threshold z» , and then above that threshold. 
The results change smoothly with z*, implying that the fi- 
nal result is not dominated by any particular cluster. Finally 
we eliminated all clusters below 6 keV (which approximately 
halves the sample in each Monte-Carlo realization). While 
the errors increase in this process, the central value only 
increases by 3-5 per cent (depending on Qm)- 

In earlier cluster samples the shape of the temperature 
function was not obviously well fit by the theoretical predic- 
tions. Much of this discrepancy has now disappeared thanks 
to better data. The question remains however as to how 
much the shape of the temperature function affects the fit, 
rather than for example the overall amplitude at some in- 
termediate temperature. To address this we have calculated 
aci for flat cosmologies by matching the observed n(> fcTx) 
for Tx = 6keV. We found values 2.5-4.5 per cent lower than 
fitting to the whole temperature function for Tx > 3.5 keV. 
Thus the shape of the temperature function is being used 
in the fit, and shifts the best fitting aci by a small, but 
non-negligible amount. 

In Fig. 7 we compare our results with those of Viana & 
Liddle (1999) and others for flat cosmologies as a function 
of fi M .§ While the Viana & Liddle (1999) best fit is within 
our error bars, the shape of the two curves is quite different. 

§ Many authors have found similar results which slightly differ 
from Viana & Liddle. We choose to focus on their study for com- 
parison, since they were very explicit about the details of their 



Figure 8. The normalization as for flat cosmologies with T = 
0.23 as derived here (solid red line) but with /3 = 1.17, compared 
with the Viana & Liddle (1999) results (dark blue dot-dashed 
line) and also compared with their approach but without the in- 
tegration over formation redshift (blue dashed line). 

This discrepancy can be traced to a number of factors: a 
newer data compilation; our likelihood fit to the entire clus- 
ter data; our use of the Sheth & Tormen (1999) universal 
mass function rather than the PS theory; not integrating 
over formation redshifts; different M—T normalization and 
scatter; and different sources of errors (in particular Viana 
& Liddle 1999 included errors in V). 

We found a relatively symmetric error bar on our final 
results for aci and <78, while Viana & Liddle (1999) had a 
very asymmetric error. It appears that this is mainly due to 
their assumed skew distribution for F. Another important 
difference is that Viana & Liddle (1999) adopted a lower 
value of (3 with respect to ours. If we chose f3 = 1.17, we 
would find the same normalization at high Qm (see Fig. 8). 
Approximately half of the discrepancy at lower Qm comes 
from the fact that we have not integrated over 'formation' 
redshift, while Viana & Liddle (1999) included such an inte- 
gration (see Fig. 8). We argued in §2 that including such an 
integration in addition to the scatter in the M—T relation 
was effectively double counting some of the scatter. 

Use of the Sheth & Tormen (1999) or Jenkins et 
al. (2001) mass function, rather than PS, lowers our results 
by 4-8 per cent (depending on Qm)- Note that the overall 
effect of the different mass function may depend on redshift, 
so that if formation redshift is integrated over then the dis- 
crepancy may even be in the other direction. 

A major source of theoretical uncertainty remains the 
value and distribution of f3 in the M—T relation of Eq. (13). 
We show in Fig. 9 how aci changes as P is scanned from 0.9 
to 1.5 (with no systematic error) in our reference model with 
Qm = 0.3 and Qa = 0.7. It is clear that a better estimate of 
the M — T relation would greatly reduce the errors on aci- 
It is also clear that to obtain an unbiased estimate of aci 

procedure. See Wang & Steinhardt (1998) for a detailed discus- 
sion of differences between some other studies. 
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Figure 9. The largely r-independent normalization u(Rci) for 
different values of (3 in the M—T relation. The error bars indicate 
the standard deviation of cr(i?ci) arising from all uncertainties 
except for the systematic error in /3, which is set to zero here. 
Our adopted central value of f3 is 1.3. 



it is essential to use an appropriate central value and error 
on (3, which is why we have included both a statistical and 
'systematic' error in our fits (see also Table 6). 

An estimate of the uncertainty introduced by various 
factors is shown in Table 6. To obtain the estimates of er- 
ror budget in this table we ran our code with the sources of 
errors described in the first column. The first case for exam- 
ple is the complete calculation, with all sources of error in- 
cluded. The second case has only the statistical uncertainty 
in the M—T relation, while the uncertainty in the value of 
the cluster temperature and the systematic scatter in the 
M—T relation were not used. The other cases show what 
happens with one or other source of uncertainty excluded. 
The mean quoted is the mean of the distribution obtained 
in each case, while the the errors are obtained by integrating 
the normalized distribution until 34 per cent of the area is 
reached on each side of the mean. 

We find that our uncertainty in erg is dominated by sys- 
tematics in the M—T relation (see also Voit 2000). The sta- 
tistical scatter in the M—T relation has an almost negligible 
effect on the central value and the error. The temperature 
errors themselves give a skewed distribution to cr 8 , but this 
effect is sub-dominant. 



6 CONCLUSIONS 

The local abundance of rich clusters of galaxies currently 
provides one of the strongest constraints on the normaliza- 
tion of the present day dark matter power spectrum on a 
scale of ~ 10 Mpc. While there have been numerous detailed 
studies of this constraint in the past, recent developments 
in both theory and observation have made it worthwhile to 
revisit this quantity. 

We have calculated the constraint on <jg arising from 
a new local sample of X-ray clusters with ASCA temper- 
atures. We have incorporated the universal mass function 



Table 6. For a flat cosmology with Qm = 0.3 we estimate how 
each source of uncertainty affects the mean and the error of our 
trrjl estimates. M— T\ Bya and M— T| s t a t represent the systematic 
and statistical errors in the M—T relation, while T represents the 
error in the temperature of the clusters. The +lcr and — lcr limits 
are found from integrating the distribution of 1000 Monte Carlos 
in each case. 



Source of error 


mean u(Rci) 


+lcr 


-1(7 


all 


0.581 


0.049 


0.050 


M-T\ sys 


0.575 


0.047 


0.047 


M-T| stat 


0.570 


0.002 


0.002 


T 


0.586 


0.005 


0.018 



determined from recent large N-body simulations, which is 
sufficiently different from the PS theory that the cosmologi- 
cal constraints inferred from cluster abundances change. We 
tried to carefully define the relation between mass and X- 
ray temperature for galaxy clusters, based on the results of 
many different hydrodynamical simulations. Another major 
difference with previous work was that we performed all the 
comparisons at the 'observed epoch' rather than carrying 
out an integration over 'formation times'. We also consid- 
ered general combinations of Um and Qa, including closed 
models. Our results are best presented in terms of a nor- 
malization at the characteristic scale for our sample, Rci, 
where the normalization is largely independent of the power 
spectrum shape. 

In the near future we anticipate a dramatic increase 
in our knowledge of the cosmological parameters from 
CMB anisotropy missions. In particular it has been fore- 
cast (Eisenstein, Hu & Tegmark 1999) that MAP will be 
able to determine erg to 14 per cent, at the same time as con- 
straining a suite of other parameters. The cluster abundance 
constraint is currently uncertain only at the ~ 10 per cent 
level, and so certainly is a pivotal limit on models - it is also 
an entirely independent constraint (at low z, in the mildly 
non-linear regime) compared with the CMB determinations 
(derived from the purely linear regime at z ~ 1000). The 
uncertainty in the derived erg could be dramatically reduced 
with improvement in the M—T relation, as well as through 
new X-ray data coming from Chandra and XMM-Newton. 

Future X-ray surveys, which go much fainter will lead 
to an increase in the size and quality of the X-ray data. We 
showed in Fig. 4 that there is lots of volume that was not 
probed at, say, 5keV. For example the ROSAT Bright Sur- 
vey (Schwope et al. 2000) contains an approximately com- 
plete flux-limited sample of 302 clusters at |6| > 30°, most of 
them lying at low redshift, which would be an excellent data- 
base if they all had good temperature estimates. More im- 
portant than just increasing the precision of the X-ray tem- 
perature measurements is an increase in the quality of data, 
for the purposes of more fully understanding the comparison 
with models. Better spatial and spectral information from 
X-ray clusters should allow more reliable estimates for the 
value of Tx which is most appropriate for comparing with 
the simulations. In a similar vein, improved optical stud- 
ies of lensing and velocity dispersions could improve mass 
determinations for individual clusters. 
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The major improvement will come through a more pre- 
cise mass-temperature relation for X-ray clusters. This will 
require both the improved data which we can foresee from 
Chandra and XMM-Newton, and improvements in mod- 
elling the gas physics relevant for understanding the X-ray 
properties of clusters. Eventually we imagine that the data 
will be fit more directly to much more ambitious simula- 
tions - this will be extremely difficult in practice, since a 
wide range of scales needs to be modelled simultaneously. 
In the meantime the piece-meal approach we have taken 
here will continue to be useful. With each of the ingredi- 
ents, particularly the mass-temperature relation, improving 
with time, the cluster abundance will continue to be a strong 
constraint on the normalization of dark matter fluctuations 
on the ~ 10 Mpc scale. 
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